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Abstract 

in : 

. We investigate multi-boson interference. A Hamiltonian is presented that treats pairs of 

bosons as a single composite boson. This Hamiltonian allows two pairs of bosons to interact 
as if they were two single composite bosons. We show that this leads to the composite bosons 
exhibiting novel interference effects such as Hong-Ou-Mandel interference. We then investi- 
gate generalizations of the formalism to the case of interference between two general composite 
bosons. Finally, we show how one can realize interference between composite bosons in the two 
atom Dicke model. 

^h' keywords: interference of bosons, quantum optics, Hong-Ou-Mandel effect. 

1 Introduction 

> 

One area in which quantum mechanics differs greatly from classical physics is with regards to 
, the statistics of identical particles. It is well known that there are two different types of quantum 

particle: bosons and fermions. The behaviour of the two types of particle is very different. A 
, characteristic feature of bosons is that they tend to bunch together. A famous example of this 

within quantum optics is the Hong-Ou-Mandel (HOM) effect [TJ[2]. This occurs when a photon 
enters each of the two inputs of a balanced (or 50-50) beam splitter. In the ideal case we find 
that the photons have a zero probability of being measured exiting both outputs of the beam 
splitter. Instead, we find both photons in one of the two outputs. One can understand this 
as being the result of destructive interference between the two possible paths that the photons 
must take if they are to be found in different outputs [3] . 

The fact that photons show HOM interference is a manifestation of the bosonic nature of 
photons. The phenomenon should, in principle, be observed in other bosonic systems. An 
interesting feature of bosons is that a system composed of several bosons is itself a boson. 
These composite bosons will also obey Bose-Einstein statistics. In principle, one should be able 
to demonstrate HOM interference with composite bosons, provided one can find an analogue of 
an optical beam splitter. 

A beam splitter is an example of an optical coupler. These are devices that allow two modes 
to interact and be transformed into two different modes. The operation of an optical coupler can 
be described by suitable Hamiltonian. For example, the Hamiltonian H = a^b+Wa can be used 
to describe a beam splitter, where a and b are the annihilation operators associated with the two 
optical modes [?]. We work throughout with a dimensionless Hamiltonian, so that our unit of 
time is defined in terms of the physical coupling constant. An important class of optical couplers 
are those that rely on nonlinear optical process. The reason for the importance of nonlinear 
optical couplers is that they have been shown to exhibit a number of novel physical phenomena 
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such as second harmonic generation [5], generation of squeezed states [S] and collapses and 
revivals of photon number oscillations [7] . We will investigate a nonlinear Hamiltonian of the 
form 

H=(tffb 2 + (tf) 2 a 2 , (1) 

where a, b are bosonic annihilation operators that satisfy [a, b] = O.We shall show that the 
Hamiltonian (fl) gives rise to novel dynamics such as HOM interference between two pairs of 
bosons. One can think of the Hamiltonian (fT} as acting like a beam splitter Hamiltonian, but 
with photons replaced by pairs of bosons. 

The outline of the paper is as follows. In section [2] we will investigate the dynamics of the 
Hamiltonian (fT}. In particular we will show that it exhibits HOM interference for composite 
bosons. The results will be interpreted in terms of a simple interference picture in section [31 We 
will then consider generalizations of the Hamiltonian (flj in section @] The case of interference 
between two composite bosons, each composed of n single bosons, is considered together with 
more general situations. A simple example of how the theory can be implemented will be given 
in section [5j Finally, we will discuss the results in section [6] 



2 A 'beam splitter' for composite bosons 

Let us investigate the dynamics induced by the Hamiltonian (fT]). It can easily be verified 
that the Hamiltonian, (fTJ), commutes with d'a + b'b, hence the total number of bosons will be 
conserved. We can, therefore, decompose the Hilbert space into invariant subspaces that have 
a fixed number of bosons. By examining the form of ([T]) we see that the state \m, n) will be 
coupled with the states |m — 2,n + 2) and |m + 2,n — 2). This implies that the subspaces 
with fixed boson number can be decomposed further into smaller invariant subspaces that have 
spanning sets of the form {\m, n), \m ± 2, n =F 2), \m ± 4, n =p 4), ...}. 

The states |1, 0), |0, 1) and |1, 1) are all eigenvectors of ([1]) that correspond to the eigenvalue 
zero. The dynamics of these states is trivial. A more interesting situation to consider is when 
we have two bosons in one mode of the field, e.g. |2, 0). It is clear that the Hamiltonian ((T|) will 
couple the states |2, 0) and |0,2). The dynamics that one observes is simply Rabi oscillations 
between the states |2, 0) and |0, 2). For a suitably chosen interaction time, we can transform the 
state |2,0) into (|0, 2) — i\2, 0))/^/2. This is reminiscent of the behaviour of a balanced beam 
splitter with respect to a single photon. The two bosons can be thought of as a single system, 
a bi-boson, which is itself a boson, albeit one with twice the energy and half the wavelength of 
a single boson [8]. 

Suppose we have a state of the form |2,2). The Hamiltonian (fTJ) will couple this state to 
the states |4,0) and |0,4). The three states |4, 0), |2, 2) and |0,4) will span a three dimensional 
subspace that is invariant under the action of equation ([I]). We shall denote this subspace as 
%3. If we diagonalize the Hamiltonian (fl} within the subspace then we find the eigenvectors 

|Ao> = i=(|4,0)-|0,4}), 

|A±) = ~(|4,0>±a/2|2,2> + |0,4)), (2) 

where H\X±) = ±4 V / 3|A±) and H\X ) = 0. The three states |4,0), |2,2) and |0,4), can be 
expressed in terms of the new basis, ([2]), as follows 

|4,0) = i(|A + ) + lA_)+V2|A >), 

|2,2) = i=(|A + )-|A_}), 

10,4} = i(|A + ) + |A_)-V2|A >). (3) 

We shall denote the evolution operator as U (t). Within the subspace ^3, the evolution operator 
has the form U(t) = exp(-iHt) = e - ,;A *|A + )(A+| + e a *|A„)(A_| + |A )(A |. 
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Suppose we have four bosons in the first mode, while the second mode is in the vacuum, i.e. 
the system is prepared in the state |4, 0). This state will evolve under the influence of (dJ to the 
state U(t)\4,0) = (V2|A ) + e~ li ^' t \\ + ) + e i4v ^*|A_))/2. We see that at time t = 7f/(8\/3), the 
state has evolved to 

U (^=) |4, 0) = \ (|4, 0) - |0, 4) - zV2|2, 2)) . (4) 

The probability of finding two bosons in each mode, at time t — ir/ (8\/3), is thus 1/2. Similarly, 
at t = 7r/(8v3) we find that the probability of finding all four bosons in the first mode is 
1/4, while the probability of finding all the bosons in the second mode is also 1/4. From the 
symmetry of the situation we see that if our state was initially |0, 4), then we would obtain the 
same probabilities for finding the bosons in the given output modes, at time at t = 7r/(8vo). 

It is instructive to compare these results to those for a balanced beam splitter with two 
photons entering one mode, while the other mode is in the vacuum. It can be shown [5] that the 
probability of obtaining a photon in each output mode is 1/2, while the probability of having 
both photons exit in a particular output mode is 1/4. The action of the Hamiltonian (pj, with 
respect to pairs of bosons, strongly resembles that of a beam splitter with respect to single 
photons. This agreement suggests that the Hamiltonian acts like a balanced beam splitter with 
respect to pairs of bosons. 

Suppose the system is initially prepared with two bosons in each mode of the field, i.e. 
the state is |2,2). We thus have one composite boson in each mode. The Hamiltonian ([I) 
will cause the state to evolve to U(t)\2,2) = (e' i4v ^*|A + ) - e iiV%t \\-)) / ^/2. The probabilities 
for detecting the bosons in a particular mode, at time t, can be calculated using the previous 
equation together with equation ([3]). A simple calculation yields the results 

P 22 (t) = |(2,2|C/(i)|2,2>| 2 =cos 2 (4V3i), (5) 
P M {t) = |<4,0|?7(£)|2,2)| 2 = |sin 2 (4V3t), (6) 

P 04 (t) = |(0,4|^)|2,2)| 2 = isin 2 (4V3t). (7) 

If we choose the interaction time to be t = n/(8y/S), then equation ([5]) implies that P22 = 0, 
while P40 = Pq4 = 1/2. The probability for detecting bosons in both output modes is thus zero. 
It can easily be shown that the output state, at t — 7r/(8v3), is just 

The term corresponding to two bosons in each output is, as expected, missing from equation (|8"j). 
This is analogous to the HOM effect, where we have a photon entering each input of a balance 
beam splitter, but where we cannot detect photons in both output modes. We should note that 
the states given in equations ^ and (|4|) are both generated using the same interaction time. 
For the fixed interaction time, t = 7r/(8v3), the action of the Hamiltonian ([1]), with respect to 
pairs of bosons, is identical to the behaviour of a balanced beam splitter with respect to single 
photons. One can thus think of equation ([I} as describing the action of a coupler that treats 
pairs of bosons as a single system, a bi-boson. Furthermore, this coupler will act as a balanced 
beam splitter with respect to bi-bosons. 

3 Interference of composite bosons 

In the previous section we showed that the nonlinear Hamiltonian (|T|) , gives rise to dynamics that 
are analogous to those observed for photons in the HOM effect. These results can be interpreted 
in terms of a simple interference picture. This interpretation will help us to understand the 
observed multi-boson dynamics. In particular, it will become clear that the zero probability of 
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Figure 1: A diagram showing the four paths the the two photons can take through the balanced 
beam splitter. 




Figure 2: A diagram showing some of the paths the the four bosons can take through the coupler. 
Part (a) shows the paths that are consistent with the two bosons being treated as a single composite 
boson. Part(b) shows one of the new paths that can occur if the two bosons can travel through the 
coupler independently. 

obtaining bosons in each mode, i.e. P22 = 0, is a manifestation of HOM interference between 
two composite bosons. 

To gain insight into the HOM effect we will first look at the familiar case where two photons 
enter a balanced beam splitter through different inputs. Each photon can be transmitted or 
reflected, hence there are four possible paths that the photons can take; this is shown in figure 
f . It can be seen that there is only one path that will lead to both photons being in a particular 
output. In contrast, there are two different paths that will result in one photon being in each 
output mode. From figure I we see that these two different paths correspond to either both 
photons being transmitted or both photons being reflected. When a beam of light is reflected 
it undergoes a phase shift of 7r/2 relative to one that was transmitted. We thus expect that 
the path corresponding to both photons being reflected will have a phase shift of 7r relative to 
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the case when both photons where transmitted. This would result in the two paths interfering 
destructively. The probability of finding a photon in each output would thus be zero, which is 
precisely what is observed in the HOM effect. 

The interference picture can be used to interpret the situation where we have bi-bosons 
entering each input of a coupler. To allow a simple comparison to the standard HOM effect we 
will refer to the bosons as photons in the following discussion. Nevertheless, the results we have 
obtained apply not only to photons but to any system of bosons. Suppose we have an optical 
coupler that has two photons entering each input. We shall assume that the optical coupler 
causes the photons to interact in the manner described by ([T]). Let us further suppose that the 
coupler treats pairs of photons (or bi-photons) as a single system. In this case we will again have 
four different possible paths for the bi-photons. Only one path will correspond to finding all the 
photons in a particular output. There are, however, two different paths that will result in two 
photons being in each output of the coupler, see figure 2 a. If we think of the bi-photons as a 
single system, then the reflection of a bi-photon will result in a phase difference of tt/2 relative 
to a bi-photon that was transmitted. The path where both bi-photons are reflected will thus 
have a phase difference of n relative to the path where both bi-photons are transmitted. These 
two paths will destructively interfere in the same manner as we described for the single photon 
case. We thus find that the probability of finding two photons in each output will be zero. We 
can see that this effect is just two boson interference, where the two bosons are composed of 
two photons. 

To emphasize that the behaviour of section [2] really is interference of two composite bosons, 
we will discuss the dynamics that occur if the two bosons are not treated as a single composite 
system. An obvious question to ask is: how many paths will result in us obtaining two bosons 
in each output of the coupler? To simplify the discussion we will again take the bosons to be 
photons. In figure 2 it is shown that there are now three paths that lead to two photons exiting 
each output- The first two paths correspond to all the photons being either transmitted or 
reflected. This is same as for the previous situation, where the photons where treated as a 
single system. The new path that we obtain occurs due the fact that the photons in each 
input can behave differently from each other. It is thus possible for one photon to be reflected 
while the other is transmitted. This will lead to an extra path that cannot be cancelled out. 
Furthermore, the two paths that correspond to all the photons being reflected and all the 
photons being transmitted, will no longer destructively interfere with each other. The reason 
for this is that each reflected photon has a tt/2 phase shift, relative to the photons that where 
transmitted. The net phase difference between the two paths is thus 2ir, which means the two 
paths will interfere constructively. We will thus obtain a non-zero probability for detecting two 
photons in each output mode. 

From the above discussion the following can be concluded. The Hamiltonian ([T]) treats pairs 
of bosons as a single system, a bi-boson. The effects we analyzed are thus examples of two boson 
interference. In particular, the fact that we cannot detect bosons in both outputs, when we have 
two bosons entering each input, is the result of HOM interference between two composite bosons. 

4 Interference of general composite bosons 

We have shown that the Hamiltonian (p} leads to HOM interference between two composite 
bosons, which are composed of two other bosons. An obvious question to ask is whether we can 
observe similar effects between composite bosons that are composed of n bosons each. The fact 
that a system composed of n bosons is also a boson suggests that we should be able to observe 
HOM interfere, provided we can find a suitable coupler. The desired action of the coupler is 
that it should treat a mode with n bosons as containing a single composite boson. The coupler 
should then act as a beam splitter with respect to these composite bosons. 

1 Actually there are four paths, however, two of these paths are identical. The indistinguishability of quantum 
particles means that these two paths should not be viewed as a separate paths. 
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Consider the following Hamiltonian 



h = (fit)»$ n + (&t)»a™, (9) 

where n is a positive integer, and a and 6 are mutually commuting annihilation operators 
that satisfy bosonic commutator relations. It can be shown that the total number of bosons 
is conserved under the action of ((9J. The Hamiltonian will thus couple the state with n 
bosons in each mode, \n,n), to the states \2n, 0) and \0,2n). These three states will thus 
span a three dimensional subspace, which is invariant under the action of (|9|). We will con- 
cern ourselves with diagonalizing (|9]) within this invariant subspace. The eigenvectors of H are 
\X' ) = (|2n,0) - |0,2n})/V2 and |A' ± ) = (|2n,0) ± V2\n,n) + |0,2n))/2, where H\X' ) = and 
H\X'±) = ±y/(2n)\\\±). A straightforward calculation shows that at time t — n/ (2^/ (2n)\), the 
probability of obtaining n bosons in each output mode, when n bosons enter each input, is zero. 
The Hamiltonian (|9|) thus gives rise to HOM interference between two composite bosons that 
are each composed of n bosons. 

Thus far we have investigated interference between two composite bosons that are each 
composed of n bosons. An interesting question to consider is whether or not one can observe 
HOM interference between two composite bosons that are composed of different numbers of 
bosons. For example, we may wish to have a field containing two photons interacting with 
another field that contains only one photon. In order to answer this question we introduce the 
following Hamiltonian 

H = {tf) n b k + (fr) k a n , n,keN:n,k>0, (10) 

where a and b are bosonic annihilation operators that act on the first and second modes respec- 
tively. Hamiltonians of the form of (flU)) have been studied before PJjJ HJJ and can be used in 
models of nonlinear couplers [T2l [13l [T4] . It can be verified that M — ka'a+nb'b commutes with 
(fTU|) . The systems Hilbert space can thus be decompose into subspaces with a fixed eigenvalue 
for M. For example, the states \n, k), |2n, 0) and |0, 2k) are eigenvectors of M that correspond 
to the eigenvalue 2nk. 

Suppose we prepare a system with n bosons in the first mode and k bosons in the second, 
i.e. it is prepared in the state \n, k). The general Hamiltonian (fT0|) will couple the state |n, k) to 
the states \2n,0) and |0, 2k). These three states will span a three dimensional subspace that is 
invariant under the action of (|10[) . We will confine ourselves to this three dimensional subspace. 
The structure of the problem is clarified somewhat if we represent (fT0|) in matrix form. We set 
|2n, 0) = (1, 0, 0) T , \n, k) = (0, 1, 0) T and |0, 2k) = (0, 0, 1) T , which yields the following form for 

mm 



H 



x 0' 
x y 
y 



(11) 



where x and y are 



n! ' V fc! 
One can easily check that the eigenvectors of (ITU) are 

|Vo) = i(y|2n,0)-x|0,2fc)), 



(2n)!fc! y {2k)M_ {u) 



|^±> = -^(x|2n,0)±A|n,fe)+i/|0,2fc)), (13) 



where A = y/ x 2 + y 2 . It can also be seen that H\if>o) = 0) H\ip±) — ±\\ip±). We can expand 
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the three vectors \2n,0), \n,k) and |0, 2k) in terms of the eigenvectors (| 1 3 1) 



|2».o) = j (y\^o) + + 

\n,k) = -j= (|tf+> - |V-}) , 

|0,2fc) = I^l^ + ^l^.)-^)). (14) 

The system will evolve from the initial state \n, k) to the state U(t)\n, k) at time t. If we 
place detectors at each output, then we will obtain one of the following three outcomes. We 
detect 2n bosons in the first mode, we detect 2k bosons in the second mode, or we detect n 
bosons in the first mode and k bosons in the second mode. The probabilities of these three 
outcomes can be easily calculated and are found to be 

2 

Pinfiit) = |(2n,0|L^)Kfc)| 2 = -^^sm 2 (Ai) (15) 

x z + y z 

Po.2k = |(0,2fc|J7(i)|n,fc)| 2 = —^sm 2 (Xt) (16) 

x z + y z 

P n ,k = \(n,k\U(t)\n,k)\ 2 = cos 2 (At). (17) 

If we choose our interaction time to be t = tt/(2X), then the probability of obtaining bosons 
in both output modes, given in equation (|17[) . will be zero. This is again an example of HOM 
interference and can be explained using the same interference picture that we described in 
section [3l One difference between the current situation and those discussed earlier is that now 
the probabilities for finding all the bosons in a given output are different for each output. Using 
equations (fT5|) and (p~6|) . we find that at t = 7r/(2A), the probabilities for finding all the bosons 
in a single output are 

P 3 „, = w » md P , 2t = JV M!=l, (18) 

n\ k\ 

where N — l/(x 2 + y 2 ). If we ignore the normalization terms in equation (|18l) . then we see 
that the probabilities are combinatorial terms that describe the number of different ways of 
distributing n and k indistinguishable particles between two modes |16) . 

We have diagonalized (fTU|) within the subspace spanned by the vectors |2n,0), |n, k) and 
|0, 2k). We could also look at different invariant subspaces. For instance the three dimensional 
subspace spanned by \2n + s, r), \n + s,k + r) and \s, 2k + r) will also be an invariant subspace 
provided r and s are integers within the ranges < s < n and < r < k. If we set \2n + s,r) = 
(1, 0, 0) T , \n + s,k + r) = (0, 1, 0) T and \s, 2k + r) = (0, 0, 1) T , then our Hamiltonian, ([T0]), will 
have the same form as (jlip . but where x and y have the values 



x 



l (2n+s)\(k + r)\ , / {2k + r)\(n + s)\ 

(n + *)!r! ^ V = V (n + r)! S ! " (19) 



The dynamics within this subspace is formally the same as in the previous case. In particular 
the probability of finding n + s bosons in the first mode and k + r in the second mode, has the 
same form as (fT7|) . A straightforward calculation shows that 

P n +s,k+r(t) = \(n + Sl k + r\U{t)\n + s,k + r)\ 2 = cos 2 {Xt). (20) 

If we choose the interaction time to be t = 7r/(2A), then we find that P n + St k+r(t) — 0. We can 
interpret this effect as an example of HOM interference between two composite bosons that are 
composed of n and A; bosons respectively. 

All of the examples that we have considered have consisted of two modes. We can, however, 
develop an analogous theory for systems with more than two modes. Wc will demonstrate this 
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be investigating a system that is composed of three modes. Suppose we have three modes that 
are coupled by a Hamiltonian of the following form 

H = (ctjW + (tf) k (P) k c n , (21) 

where a, b and c are bosonic annihilation operators that act on the first, second and third system 
respectively. This Hamiltonian will couple the states \n,k,k), |2n,0, 0) and |0, 2k, 2k). If we 
set |2n,0,0) = (1,0, 0) T , |0,2fc,2fc) = (0,1, 0) T and \0,2k,2k) = (0,0, 1) T , then will be a 
matrix with the same form as but where x and y have the form 

(22) 

The dynamics of this system will thus be formally the same as for the general Hamiltonian (TIT)]) , 
within the three dimensional subspace. In particular we can observe HOM interference, i.e. we 
can choose the interaction time so that a state \n, k, k) will be transformed to a superposition 
of the states |2n,0,0) and |0,2fc,2fc). 



5 Example of Hong-Ou-Mandel interference in the Dicke 
model 

In the previous sections we have shown how assemblies of bosons will interact as a single com- 
posite system, under the action of a suitable Hamiltonian. We have not, however, discussed 
how one could realize the theory in an experimentally feasible context. In this section we will 
give a concrete example. 

Suppose we have two atoms in a cavity that are interacting with a single mode of the 
electromagnetic field within the cavity. For simplicity we shall assume that the cavity mode is 
resonant with a particular transition between two energy levels of the atoms and consider only 
this transition for each atom. The atoms can be thought of as having only two energy levels, 
the ground state, \g) and the excited state |e). This is an example of a two atom Dicke model 
|17) . We will consider only the symmetric states of the two atoms. We can thus represent the 
atomic state in the form |m), where m equals the number of atoms in the excited state. For 
instance, |0) = \g,g), |1) = (\e,g) + \g,e))/y/2 and |2) = |e,e). We can define the individual 

atomic dipole operators as follows: atf) will raise the j-th atom from the ground state to the 

(i) 

excited state, and will lower the j-th atom from the excited state to the ground state. We 

also define = 2a+&_ — 1. The total atomic dipole operators, S z , S± will be taken to be 

the sum of the individual atomic dipole operators, hence S z ,± = p^Iz + "z~±- One can verify 
that 

S + \m) = y/2\m+i), 
S-\m+l) = V2|m), (23) 

where m = 0, 1. 

A useful way of describing this system is to make use of the Schwinger representation for 
angular momentum |18j . In this approach we introduce two fictitious bosonic modes b\ and b 2 . 
We then make the following substitution for the atomic dipole operators 

S+ = bxb\, 
S- = b\b 2 

S z = (bib 2 -b\k)/2. (24) 

The atomic states can be represented by Fock states \j,f), where |m) = |2 — m,m). The 
Fock state |1, 1) is thus equivalent to the symmetric state (\e,g) + \g,e))/y/2. The Hamiltonian 
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that describes the interaction between the two atoms and the field mode will have the form 

[322 EES [13 

Hi = ahbl + a r b\b 2 , (25) 

where a is the annihilation operator for the cavity mode. We can see that the Hamiltonian (1251) 
has the same form as (f2"Tj) . where n = k = 1. The total state of the cavity field and the atoms 
can be expressed in the form \N,j,j'), where iV is the number of photons in the cavity mode and 
\j,j') is the Fock state that represents the state of the atom in the Schwinger representation. 

The cavity field will be prepared so that it contains a single photon, while the atom pair 
is initially prepared in the symmetric state |1, 1). The initial state of the full system is thus 
1 1 , 1 , 1) . After the atoms and cavity mode have interacted for a time t = 7r/(2A), the state of 
the system will have been transformed to 




The state (l2"6l) does not contain the term |1, 1, 1). We can view this as HOM interference between 
two composite bosons, where the first composite boson consists of a photon and an atom, while 
the second boson is composed of a single atom. 

6 Conclusion 

We have investigated a special type of multi-boson interference. A Hamiltonian, ([1]), was pre- 
sented that treats two bosons as a single composite boson. We have showed that the Hamiltonian 
allows two composite bi-bosons to display HOM interference. The interpretation of this effect 
was corroborated by a simple interference picture. In particular, we showed that the observed 
phenomenon was a direct consequence of the two bosons acting as a single composite boson. If 
the two bosons had acted independently, then we would have observed bosons in both outputs. 

The task of engineering a coupler that causes the bosons to interact in the manner described 
by ([TJ is non-trivial. In practical situations we could have higher order terms appearing in 
the Hamiltonian. One might expect that these terms will always lead to the interference being 
washed out. This, however, is not always true. As an example, consider adding the higher order 
terms Afa 4 ^) 4 + (a^) 4 6 4 ] to the Hamiltonian (p}. In this situation the states |4, 0), |2,2) and 
1 , 4) are still coupled. Furthermore, it can be shown that we still observe HOM interference, 
but where now the interaction time needed to observe this will have changed. 

We have shown that within the 2 and 4 boson subspaces, the Hamiltonian (JT]), acts like a 
beam splitter with respect to bi-bosons. This will not hold within higher dimensional subspaces. 
It can be shown that the dynamics of N bi-bosons, N > 2, is not the same as the dynamics of 
N photons interacting via a beam splitter Hamiltonian. This is because the matrix elements of 
two Hamiltonians are different within the corresponding invariant subspaces. 

An assembly of N bosons can be viewed as a single composite system, which is also a 
boson. This observation means that it is possible to observe two boson interference between two 
general composite bosons, provided one can find a suitable interaction. We presented a general 
Hamiltonian, (fTU|) . that would allow one to observe interference between composite bosons that 
where composed of n and k bosons. We were able to observe HOM interference between these 
two composite bosons. This general type of two boson interference exhibited interesting effects 
such as the following. For n = fc, the probability of obtaining all the bosons in the first output 
was the same as the probability of obtaining all the bosons in the second output. When n =/= k, 
we find that these two outcomes occur with different probabilities. The reason for this was 
that we observe different numbers of bosons in each output, either 2n or 2k. The probabilities 
for observing 2n or 2k bosons in a given output are, up to a normalization constant, just 
combinatorial terms that give the number of ways of removing n or k indistinguishable bosons 
and then redistributing k or n new bosons with the remaining bosons in the given output. 

Finally, we showed how the two atom Dicke model could give rise to dynamics that could be 
interpreted as HOM interference between two composite bosons. In this case the interference was 
between an atom and photon pair and a single atom. In principle one could observe interference 
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between two composite bosons in other systems, such as Bose-Einstein condensates, the N atom 
Dicke model or between phonons in trapped ions. The implementation in terms of phonons is 
of particular interest as the vibrational dynamics of trapped ions operating within the resolved- 
sideband regime, can be described by a nonlinear Hamiltonian [2U [22] ■ 
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